REPORT No. 121 



THE MINIMUM INDUCED DRAG OF AEROFOILS 

By MAX M. MUNK 
National Advisory Committee for Aeronautics 



8077 J-Jl 1 




111. />£ : ^ 



' 1 '«■ 'J 




i 



REPORT No. 121. 



THE MINIMUM INDUCED DRAG OF AEROFOILS. 

By Max M. Munk. 



, V , INTRODUCTION. 

The following paper is a dissertation originally presented by the author to the University 
of Goettingen. It was intended principally for the use of mathematicians and physicists. The 
author is pleased to note that the paper has aroused interest in other circles, to the end that 
the National Advisory Committee for Aeronautics will make it available to a larger circle in 
America. The following introduction has been added in order to first acquaint the reader 
with the essence of the paper. 

In the following development all results are obtained by integrating some simple expressions 
or relations. For our purposes it is sufficient, indeed, to prove the results for a pair of small 
elements. The qualities dealt with are integrable, since, under the assumptions we are allowed 
to make, they can not be affected by integrating. We have to consider only the relations 
between any two lifting elements and to add the effects. That is to say, in the process of inte- 
grating each element occurs twice — first, as an element producing an effect, and, second, as an 
element experiencing an effect. In consequence of this the symbols expressing the integration 
look somewhat confusing, and they require so much space in the mathematical expression that 
they are apt to divert the reader's attention from their real meaning. We have to proceed up 
to three dimensional problems. Each element has to be denoted twice (by a Latin letter and 
by a Greek letter), occurring twice in a different connection. The integral, therefore, is sixfold, 
aix symbols of integration standing together and, accordingly, six differentials (always the same) 
standing at the end of the expression, requiring almost the fourth part of the line. The mea nin g 
of this voluminous group of symbols, however, is not more complicated and not less elementary 
than a single integral or even than a simple addition. 

In section 1 we consider one aerofoil shaped like a straight line and ask how all lifting 
elements, which we assume to be of equal intensity, must be arranged on this line in order to 
offer the least drag. 

If the distribution is the best one, the drag can not be decreased or increased by transferring 
one lifting element from its old position (a) to some new position (&). For then either the 
resulting distribution would be improved by this transfer, and therefore was not best before, or 
the transfer of ah element from (b) to (a) would have this effect. Now, the share of one element 
m the drag is composed of two parts. It takes share in producing a downwash in the neighbor- 
hood of the other lifting elements and, in consequence, a change in their drag. It has itself a 
drag, being situated in the downwash produced by the other elements. 



:€taiidflriiig only two elements, Fig. t shows that in the case of the lifting straight line the 
two downwashes, each produced by one element in the neighborhood of the other, are equal. 
For this reason the two drags of the two elements each produced by the other are equal, too, 
and hence the two parte of the entire drag <# the wings due to one element. The entire drag 
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produced by one element has twice the value as the drag of that element resulting from the 
downwash in its environs. Hence, the entire drag due to one element is unchanged when the 
element is transferred from one situation to a new one of the same downwash, and the distribu- 
tion is the best only if the downwash is constant over the whole wing. 

In sections 2 to 6 it is shown that the two parts of the drag change by the same value in 
all other cases, too. If the element* Are situated in the game transverse plane, the two parts are 
equal. A glance at Fig. 2 shows that the downwash produced by (1) at (2), (3), (4), and (5) 

: '■ ' - # 

is equaL But then it also equals the downwash due to (4), say, produced at (1). This holds 
true even :f or the component of the downwash in the direction of the lift if the elements are nor- 
mal to each other (Fig* 3.) ; for, this component is proportional ^y/r*, awarding to the symbols 

.:>:,. \ \r ?v* . *<t , >' ^ , - Jj ::..i. : > • ru 
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of the figure* , Hence, it is proved for lift of any inclination, horizontal and vertical elements 
being able, by combination, to produce lift in any direction., \ ; , 

There remains only the question whether the two parts of the drag are also equal if the 
elements are situated one behind the other — that is to say, in different longitudinal positions. 
They : are not; tut their sum is independent of the longitudinal distance apart. To prove 
tiuS| add in liig. 4 to the lifting element (2) a second inverse lifting element (3) with inverae 

. 7 • • • ' • i, . l-» <-.••.• •• r ': ' 

linear longitudinal vortices in the inverse direction. The reader observes that the transverse 
vortices (2) and (3) neutralize each other; the longitudinal linear vortices, however, have the 
same sign, and all four vortiees f onn a pair of vortices running from infinity to infinity. The 
drag, produced, by the combination of (1) . and this pair, is obviously independent of the longi- 
tudinal positions of (1) and (2). But the added element (3) has not changed the drag, for (1) 
and (3) are situated symmetrically and produce the same mutual downwash. The direction 
of 4he lift, however, is inverse, and therefore the two drags have the inverse sign, and their sum 

is zero. . . . ' , . . • ' * ;'. , '7 '/< * 

If the two lifting elements are perpendicular to each other (chapter 5), a similar proof can 
be given. 

Sections 6 and 7 contain the conclusions. The condition for a minimum drag does not 
depend upon the longitudinal coordinates, and in order to obtain it the downwash must be 
assumed to be constant at all points in a transverse plane of a corresponding system of 'aero- 
foils. This is not surprising; the wings -act like two dimensional objects accelerating the air 
passing in as* infinite transverse plane at a particular moment. Therefore the calculation 
leads to the consideration of the' two dimensional flow about the projection of the wings on a 
transverse plane. 
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Section 8 gives the connection between the theory in perfect fluids and the phenomenon 
in true air. It is this connection that allows the application of the results to practical questions. 

1. THE LIFTING STRAIGHT LINE. 

A system of aerofoils moving in an incompressible and frictionless fluid has a drag (in the 
direction of its motion) if there is any lift, (perpendicular to the direction of its motion). The 
magnitude of this drag depends upon the distribution of the lift over the surface of the aerofoils. 
Although the dimensions of the given system of aerofoils may remain unchanged, the distribu- 
tion of the lift can be radically altered by changes in details, such as the aerofoil section or the 
angle of attack. The purpose of the investigation which is given in the following pages is to 
determine (a) the distribution of lift which produces the least drag, and (6) the magnitude of this 
toinitnum drag. ' ' 

Let us first consider a single aerofoil of such dimensions that it may be referred to witn 
sufficient exactness as a lifting straight line, which is at right angles to the direction of its flight. 
The length or span of this line may be denoted by I Let the line coincide with the horizontal, 
or x axis of a rectangular system of coordinates having its origin at the center of the aerofoil. 
The density of the lift 

■ » 

where A, the entire lift from the left end of the wing up to the point x, is generally a function of 
* and may be denoted by/ (x). Let the velocity of flight be v„. 

The modern theory of flight 1 allows the entire drag to be expressed as a definite double 
integral, if certain simplifying assumptions are made. In order to find this integral, it is neces- 
sary to determine the intensity of the longitudinal vortices which run. from any lifting element 
to infinity in a direction opposite to the direction of flight. These vortices are generally 
distributed continuously along the whole aerofoil, and their intensity per unit length of the 
aerofoil is r * 

r __L . if (2) 

where p is the density of the fluid. Now, for each lifting element dx, we shall calculate the down- 
wash w, which, in accordance with the law of Biot-Savart, is produced at it by all the longi- 
tudinal vortices. A single vortex, beginning at the point x, produces at the point x-{ the 
downwash 

dw „ 1 . dA ' . (3) 
4irpv 0 {-X 

Therefore the entire downwash at the point {is 

ft* 



(4) 



The integration is to be performed along the aerofoil; and the principal value of the integral m 
to be taken at the point x - {. This rule also -applies to all of the following integral*. Hence * 
follows that the drag according to the equation ''" 

• A' (5) 

ax . v% 



I 

"1 



(6') 



, 5^ U prendtl, Tr»gfltt«etth*orie, I. Mitteilung. N&chrichUn der G*. <L Wiaa. «i Gflttingen, 1918. 
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or, otherwise expressed,' < f ■ • • " : ■rj[-....- >. t . -V^ : .-.\;-. ■ ,-'., v: .< •> 



/* here fignifies the derivative of / with respect to x or {. The entire lift is represented by 



V"'-- ; '<7> 



Hence the solution of the problem to determine the best distribution of lift depends upon 
the determination of the function/ so that the double integral ; v . 



shall have a value as small as possible; while at the same time the value of the simple integral 

•^i- fj (x)<fe-oonst * . a (9) 

k fixed. r * -: J • ,• ■ v7 •, - >; ^ •. 

The first step towards the solution of this problem k to form the first variation of V t 

, The second integral on the right side of (10) can be reduced to the first. By exchanging the 

symbols x and { and by partial integration with respect to x, considering/*® as the integrable 
factor, there is obtained 

The second member disappears since /-<? at the limits of integration* Further, the right 
hand part of (11) 



upon substitution of the new variables x and t-z~i for jj and f, is transformed into ' 



5 

; «H 

»lf thk wtn not tnw, thai* would bt taSaiM rdcdtiM tt Omm potato. 
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Now- 



or, since /disappears at the limits of integration, 



**4 s+ i 



which, upon the replacement of the original variables, becomes 



CfiQ 



so that, finally, 



J/^J^}-f{a/(x)irfe*} (12) 
-4 -r - ~* 

Substituting this in (10) there finally results < < • ' 



From which the condition for the minimum amount of drag, taking into consideration the 
second condition (9) , is 



4 



or, when equation (4) is taken into consideration . 

, ': w^eonsfc-ttfc - ( 15 ) 

The necessary condition for the minimum of drag for a lifting straight line is that the down- 
wash produced by the longitudinal vortices be constant along the entire line. 

That this necessary consideration is also sufficient resulte from the obvious meaning of 
the second variation, which represents the infinitesimal drag produced by the variation of the 
lift if it alone is acting, and therefore it is always greater than zero. 

r . 2. PARALLEL LUTING ELEMENTS LYING IN A TRANSVERSE PLANE. 
- ' The method just developed may be applied at once to problems of a more general nature. 
If, instead of a single aerofoil, there are several aerofoils in the same straight lme perpen- 
dicular to the direction of flight, only the limits of integration are changed m ^^^P^- 
The integration in such cases is to be performed along all of the aerofoils. However, this is 
none^ential for all of the equations and therefore the condition for the minimum drag (equa- 
tion 15) applies to this entire system of aerofoils. : 
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Let us now discard the condition that all of the lifting lines are lying in the same straight 
line, but retain, however, the condition that they are parallel to each other, perpendicular to 
the line of flight as before, and that they Are all lying in a plane perpendicular to the line of 
flight. Let the height of any lifting! line be designated by z or f . Equation (3) transforms into 
a similar one which gives the downwash .produced at the point z } 2 by the longitudinal vortex 
beginning on the lifting element at the point f : 

The expression, which must now be a mfrmn imj is \ | ^ 

with the unchanged secondary condition 

•^j-jyt** *)^- const. ■ - (9a) 

These integrals are to be taken over all of the aerofoils. 
This new problem may be treated in the same manner as the first. 



is always to be substituted for — • It may be shttwtt that this wibetitulian does ^«>l 
affect the correctness of equations (10) to (15). Therefore 



w~const-tfl, (15a) 

is again obtained as the necessary condition for the minimum of the entire drag. 
* Finally, this also holds true for the limiting case in which, over a limited portion of the 
transverse plane, the individual aerofoils, like Venetian blinds, lie so closely together nthat 
they may be considered as a continuous lifting part of a plane. Including all cases which 
have been considered so far, the condition for a mmm>iim Q f fa^g can \> t stated: 

Let the dimensions of a system qf aerofoils be given, those in (he direction of flight being small 
in comparison with those in other directions. Let the lift be everywhere directed vertically. Under 
these condition*, the downwash produced by the longitudinal vortices must be uniform at all points 
on the aerofoils in order that there may be a minimum of drag for a given toUdiift^, y < v 

' < j «. THREE DIMENSIONAL PARALLEL LIFTING ELEMENTS. 

The three-dimensional problem may be based upon the two-dimensional one. Let now 
the dimensions in the direction of flight be considerable and let the lifting elements be dis- 
tributed in space in any manner. Let y or 1 be the coordinates pf any point in the direc- 
tion of flight. For the time being, all lifting forces are assumed to be vertical , 

The calculation of the density of drag for this case is somewhat more complicated than i* 
the preceding cases. Consideration must be given not only to the longitudinal vortices, which 
are treated as before, but also to the transverse vortices which run perpendicular to the lift at 
any pqint and to the direction of flight. Their intensity at any point where there is a lifting 
dementis «,v < 



VP t^P ' 

The density of drag, W 1 now has two components, % and W t , the first heing due to the trans- 
verse vortices and the second to the longitudinal vortices. ■■■■■.■'>< ■ 
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For the solution of the present problem only the total drag of all lifting elements 

W-fWdx 

1 is to be considered. In the first place it will be shown that the integral of those parts of the 
density resulting from the transverse vortices 

Wi-fW^dx 

does not contribute to the total drag. A small element of one transverse vortex of the length dx 
at the point (x, y, z) produces at the point ({, f ) the downwash 

^-^^/(x, y,z)-dz (16) 

where 

r»- (*-*)•+ Oi-lrt'+tf-*)'. 

Therefore 

lfi-J5^JJ/&.r»t)./tt f f,f) ^dzdi. (17) 

This integration is to be extended over all the aerofoils. It is possible to write this expression 
in such a manner that it holds for a continuous distribution of lift over parts of surfaces or in 
space. This is true, moreover, for most of the expressions in this paper. Now, exchanging the 
variables x, y, z, for £, if, f , in equation (17) does not change the value of the integral, since the 
symbols for the variables have no influence on the value of a definite integral. On the other 
hand, the factor (* - y), and therefore the integral also, changes its sign. Hence 

W;--P a -0 (18) 

and, as stated, 

F- W % . (19) 

Therefore the entire drag may he calculated without taking into consideration the transverse 
vortices. 

, The method of calculating the effect of the longitudinal vortices can be greatly simplified. 
At the point (£, if, f ) that part of the density of drag resulting from a longitudinal vortex begin- 
ning at the point (x, y, z) is 

Wi'-j^/tt. *3* (20) 

where 
and 

The entire drag is _ ; ; , 

'*-f**'&-£jJf<if n , t)f (x, y, z) tftftfx. (22) 

Now, in the double integral (22) the variables x, y 9 z may be exchanged with iy, f , as before, 
without affecting the value of the definite integral. Partial integration may then be performed 
twice, first with respect to { and then with respect to x. The substitution results in 



w ~£~JS f ( *' y > z)f ({ ' * f) ***** (23) 



S077*-21- 
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fis obtwke^fc * upon the exchange of variables. Its "value is therefore 

- #-Af^&>-V-*> , +<»-y> , +<^-*> , •• ( 

When partially integrating with respect to dt, the integrahle factor is/ ft, v, t) 



In the subsequent partial integration with respect to dx, the integrahle factor is ^ -gr 
Finally, by addition of (22) and (26), there is obtained 



»Vf ■■■.!! t 



v tPj j . , • , , 



tin;-?" .10/- 



"i may now be substituted in (24) for the Variable of integration ».'' Then t changes to 

I and with the exception of the sign the integrand in (21) agrees with the resulting one in (22) 

•>.l- -. ••• • ■ ;•.:!•=< " ■■ >, '• . " ' 

Subtracting (28) from (21) there results finally ^ ; „• . .>i;»,j* 

«>>i> , . • ir i» 

;iV«ii:i*iin-- -•iji ; --,r" •*•» n,v> > ! ii."i<r f'ti'ihuti. >.- s "<?J 1-* J»vlb ■ i^iiJiJi^'f:^ T«i.hr.vi'n'n i '»i!T 
- Bfanw, and. therefore the entire right side, qf, equation ; JSiieen, to.be ^ mdj^ent pf the 
longitudinal coordinates y of the lifting elements. , : f- . .' ;H.<; ' ' " 

Therefore the entire resistance of a three-dimensional system of aerofoUs imth paralUl Upmg 
elements does not depend upon the, longitudinal position? of the lifting elements. 

4. LIFTING ELEMENTS ARRANGED IN ANY DIRECTIONS IN A TRANSVERSE PLANE.,,..! 

The problem considered in section 2 can also Joe generalized in another way. For the present 
the condition that all lifting elements be in one transverse plane may remain. However, they 
need no longer be parallel, and the lift may be due to not only a great number of infinitesimal 
lifts dA but also to similar transverse forces dB. In the first place let the direction of all lifting 
dements be arbitrary, but such that there is* minimum d>i»g, and let this direction be an 
unknown quantity to be determined. '__•„...> . •••• . 

In the present problem it is desirable to consider a continuous distribution Of IrftoVdr given 
Areas instead of lines. , The, last case can be deduced from the first at any time by passing to 
the limit. . . , . . , i U . .. " ' ' ' 



Let A'=f(x,z) be the density of the vertical lift per unit area, and B' - F(x, z) the density 
! of tne lkteral force per unit area. The lateral force is coiwiderod posil3ve %heii' acting m the 
''positive direction of the ' X-axis. Then the density of the ;trattsverse Vorfices has the 
nonents — *A' and — - -B'. The density of the longitudinal vortex is the divergence of the 

^density of the tommnil^ «r ^(^-^). i.Tte'"' Wtudii^_ Torti^ _ beginning 



minimum induced dba& o* aebofoil& m 
at the point ir, z) therefore produce at the point ($;?) the downwash and the transverse velocity 

^-sbGsrf)**^ - : '<»> 

According to the above, the density of the drag is ' t: " L H ' " ; 1 *' 1 1 • ^ r; J - ' 

With these symbols there results for the tolal drag the expression 



• .i -v - - - ' ' 



All of these integrals are to be taken over all of the lifting surfaces. Now the first two 
integrals have forms corresponding to the integral in (S), and therefore there is a possibility of 
substituting for these. . A similar relation also holds for the last two integrals. . . For '^Am- 
ple, the variation of the third integral is ,^v. ti > 



Now in the first term on the right-hand side, the variables 4 an4 z may bf exchanged with £ 
: and It may then be partially integrated with respeet to*<?$, the integrable factor being df 1 (£,f ) . 
Ogives . , - * \ • . ... , J j 'i Vji 



This may be partially integrated with respect to the integrable factor being 

Hence the first term of the variation of the thitd mt^gial.df b^) can be transformed into the 
second term of the variation of the fourth integral of this equation. In a similar, manner t the 
two other terms may be transformed into each other. It js therefore demonstrated that the 
variation of the entire drag may be written. £5 > , V s ^ * | *, , \ 

iW^2jjSf^dzdz+2fjBF'U : dxd2 (13b) 
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Two problems of variation can now be stated. In the first plaoe limited parte of the surf aces 
may be at our disposal, oyer which the vertical lift A and the horizontal transversal force B 
may have any distribution. Only the total lift 

A-Jjf(zf) daxfe-const. (9b) 
will be given in this case. f 

. Then . • " i\x : : \ : - ■";'-}// 

to •* const. ~ tt— o (15b) 
is the condition f or the least drag. ^ ; - if ^: - ' 

If, however, the lifting parts are similar to lines, there is generally one other condition to 
fulfill. It is then required that the lift disappear everywhere along the direction of the aerofoils. 
That is to say, 

/sin0-FcoB0-0 , ; (34) 

where fi is the angle of inclination of the aerofoil to the horizontal X-axis. In order to add the 
new requirement (34) a second Lagrange constant m is introduced. The condition for the least 
drag is now 

• +k+ 5h-V-aftir -. '.. ; v (34ft) 

and after the elimination of n 

w cos P+u sin p—w, cos 0 (15c) 

the constant 2 X being replaced by - u>„ as before. In words: ; v ' 

If aU lifting dements -are in me transverse plane, (he component of the velocity perpendicular 
to ike wings, produced by the longitudinal vortices, must be proportional, at aU lifting dements, to 
the cosine of (he angle of lateral indination. ? ! ,s - r 

*. LIFT DISTRIBUTED AND DIRECTED IN ANY MANNER* 

The results obtained previously can be generalized not only for lifting elements distributed 
in a transverse plane but also for lifting elements distributed in any manner in space. That 
part of the total drag resulting from the transverse vortices is, in the general case 



(17a) 



Both terms have the same form as the integral in (17). The demonstration for (17) therefore 
applies to both. In the general case also the total drag can be calculated from the longitudinal 
vortices without taking into consideration the transverse vortices: 



(22a) 



In this aa in (20), •»•>>, 
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The first two terms in (22a) have the same form as the right-hand side of (22), and the same 
conclusions are therefore valid for each. It can be proved directly for (22a) as for (22) that 
each of the two double integrals is independent of the longitudinal coordinates of the lifting 
elements. This proof can now be extended over the last two integrals of equation (22a), 
The third integral, after changing the variables, becomes 

/////// (z, y> z) F' (g, v, ft* dxdydzd&vdt- (35) 
J//////U, *, f) F' (*, y, z) * x dtfndtdxdydz 

where 

Now, let F' be chosen as the integrable factor and be partially integrated with respect to z. 
' ///////<*' V, *> (fe f) (36) 

/////// fo* 0 F (x ' y- 2) s *i*wp*«M»- 

As in the previous cases, the second integral to be expected vanishes since / as well as F 

disappear at the limits of the integration. Next Jr - - ^ ^ischosen asthe integrable factor 

and partially integrated with respect to x. By by analogy, is meant 

k .■ \ . '., . -.. J 1 . *..; ( V; - /. ^ \ » . —>•*.» . • 1- f f— g j ' ■■■ 

*» — i?J 

/////// (x, y>2 ) < 37 > 

Now + a may be transformed, the variable a in the defining equation being replaced by ij + y-s. 
The result is that 

£ J *Z* e - (£- x)» + (n - «)' + (f - a) 1 - 

It is seen that the integrand agrees with that of the defining integral Therefore, and since 
the right-hand side of (37) contains the same function under the double integral as the fourth 
turn in (22a), this fourth term can be combined with the transformed third member. This 
gives ■~..Ar/:«, : -:\ ; - > ' • . . f .- : ■• 

/////// 0»,y, *) F'«,.n,f)*,4«i 1 «i^f+ (38) 

//////F (*, y, x)J'% n, ' ri-t&fyibi&rt- 

where ■: 

and therefore the two sides of (38) are independent of y. This is therefore demon- 
strated for the whole right-hand side of (22a). 
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In general it can therefore be said: - ur m:r". . .n*ff \x,v. vi ■'■■*"» vfT 

The total resistance is always independent of the longitudinal coordinates of the lifting dements; 
And further: ^ t ■ k'/y- ..; ^ju- -A -v.v . r • 

The most favorable distribution of tike lift, with reference to the total drag, occurs when this is 
also the case for the projection of tfu lifting dematie on a transverse plar^ j v m-i •; 

That is to say, all of the lifting elements are projected on a plane perpendicular to the 
direction of flight, and any element so obtained has a lift equal to the sum of the lifts of all lifting 
elements projected onto iLu*W;v'A», v ; ' * .; > a ■ i\Y*T\'\ 

6. DETERMINATION OP THE SOLUTIONS. 

The previous demonstrations show that the investigation for the distribution of lift which 
causes the least drag is reduced to the solution of the problem for systems of aerofoils which are 
situated in a plane perpendicular to the direction of flight. In addition, the condition for least 
drag (15c), which becomes the condition of uniform downwash (15) if the lift is vertical, leads 
to a problem which has often, been investigated in. the theory of two-dimensional flow with * 
logarithmic potential. The flow produced within the lifting transverse plane by the longitudinal 
vortices originating in it is r indeed, of this type. Each such vortex produces a distribution of 
velocity such as is produced by a two-dimensional vortex of half its intensity, and the whole 
distribution of velocity is obtained by adding the distributions produced by the longitudinal 
vortices. The potential flow sought is determined by the condition of (15c). Let it be com- 
bined with the flow of constant vertical upward motion — « r The resulting flow satisfies 
the condition at the boundaries > . 

m « ■ « cos sin 0-0 (39) 

and there results, for the case of lifting lines: 

The two dimensional potential flow is cf ike type that encircles ihe lifting lines, and at a great 
distance the velocity is directed upwards and has\jh4 value w^ - «v 

Within lifting surfaces the velocity is zero - according to the condition (15b), and the fluid 
therefore flows around tire contend 

The intensity of the longitudinal vortices at any point is twice the rotation of the two 
dimensional flow. In the cise of the lifting lines/ therefore, the density of the longitudinal 
vortices is double the discontinuity of velocity from one side to the other. The intensity of the 
transversal vortices is determined by integrating the longitudinal vortices along the aerofoils 
and therefore equals twice the difference of the velocity-integral produced on the two sides of 
the aerofoil. Now the integral of the velocity produced is identical with the potential and 
hence it appears: ' r * : - - u ; | ^ 

The density of the lift perpendicular to the lifting line is proportional to the discontinuity of 
potential f%—v x * and has the value ^ 

; " ; V4"+*"-2t^(*-*> " : < * * (40) 

Hence the total lift obtained by integrating over all aerofoils is 

A^v^fa-ridx ' •/ (41) 

Sometimes a transformation of this equation is useful. In order to obtain it, suppose that 
all of the lifting lines are divided into small parte. Then, on the two ends of each lifting element 
there begin two inverse longitudinal vortices, the effect of which on a distant point is that of 
double vortex. Their velocity-potential v and their stream function ^ may be combined in the 
complex function ^ + i*p, and, not considering the Existence of a parallel flow, which is without 
any importance in the calculation, this complex function has the form for a lifting line, 

r ^ *t 

i": : X ■- ; «... 
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, +V _fM±^? (42a) 



where z represents x+iy and 3 0 «x*+{y 0 r x 0 and y 0 being the coordinates of the lifting elements of 
the line. For a lift distributed over areas a similar equation can be formed. The integration of 
(42) gives 

Now the residuum of the integrand at infinity is dA + idB and therefore the residuum of the 
integral ia A + iB. Therefore the expression can be written. ; . 

■-r :> ; V f^rviV ; -0.) : / r A-2v&R [Resty +%<(>)] ; (41a) 

where the last part means the real part of the residuum of ^ + at infinity. In the most im- 
portant case of horizontal aerofoils the residuum itself is real and can be used directly to calcu- 
late the lift. The density of drag at any point is proportional to the perpendical component of 

the density of lift and isW-^-A', from which results F-^°- A Making use of (41) one obtains 



; ^ ^'^fitRes7*-ftV)3 (43a) 

The integral in the denominator of (43) represents an area characteristic of the system of 
aerofoils investigated. Frequently the easiest method of calculation is to assume from the 
beginning the velocity w 9 at infinity to be unity. 

The case of the lift continuously distributed over single parts of areas is derived from the 
preceding one by passing to the limit. Since the vertical velocity w disappears at all points in 
the lifting surfaces, the velocity is zero at all points and the rotation vanishes. 

Therefore, in the case of the most favorable distribution of lift, att of the longitudinal vortices 
from the continuously lifting areas begin at the boundaries of the areas. 

Equations (43) and (43a) remain. The distribution of lift is indeterminate to a certain 
extent. On the other hand, it is possible to connect the points of the contour having the same 
potential <p by strips of any form, and it is only necessary that the lift be always perpendicular 
to the strip and its density have a constant value along the whole strip. According to equation 
(40) this equals the difference of the potential at the contour between the two borders of the 
strip. Worthy of note is the special case in which all of the strips mm along the contour, thus 
coming again to the case of lifting lines. It appears that : 

Closed lines have the same mimimum of drag as the enclosed areas when continuously loaded. 

Especially important are those symmetrical contours which are cut by horizontal lines in 
only two points. With such the limitation to vertical lift does not involve an increase of the 
minimum drag. For this case it appears that: 

The density of the vertical lift per unit area must be proportional to the vertical component of 
the velocity of the two-dimensional flow at the point of the contour of the same height 2. It is 



The corresponding density of drag is 



d^-*,** 1 ** (45) 
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- v ^i3,r ,. : *>.>.. .. ...... j : 7. EXAMPLES OF CALCULATIONS. ■ ;,^ r 'vV.: - v,,..- 

Examples of calculation of the previous demonstrations can be based on any calculated 
two-dimensional potential flow around parts of lines or areas. The simplest flow of the first 
kind is that around a single horizontal line. It leads to the problem investigated at the begin- 
ning of this paper. 

In this case the potential is the real part of Vj^-i, where f denotes <r The lifting 
line joins the two points z-0 9 x--l And z-0, and has the leqgth 2. Th e veloc ity a* 

infinity is «?— 1. The discontinuity of potential along the lifting line is ? a — ? t -2^1 — A The 
density of lift is distributed according to the same law, therefore if plotted over the span the 
density of lift would be represented by the half of an ellipse* 

lie minimum drag is v ![ 

If, instead of the value 2, the span had Uie general value J, the minimum drag would be 

This same result has been obtained by Prof. Prandtl by another method.* 
The simplest example for a lifting vertical area is the circle. Let its center coincide with 
the origin of the system of coordinates. Then the potential of the flow around this circle is 

- ■ ■ ;■ ■■ ; ; r\- ' • (48) 

where r-V??^ At infinity Under &e Widition of and according to equation (40) 

the density of lift is . .. ri . Vj; !" ^ ^v^i^^t^'Vix ..< v: , - fi ;- 

^ i ^-2 v i ft^ii:^'^^;: ^(^ ,^ ; tar -vi; >.-., V,;, ; 

This results in a constant density of lift of A 9 -2. Therefore the drag is 



The double integral is to be taken over the circle. If the general case for the diameter 
equal to D be considered, then the least drag ia, ; ,, :M tv v.;-- - t * j 



Hence in respect to the minimum drag the circle is equivalent to a lifting line having 
a length i/2 times the diameter. " : , ' ; 

A lifting circular line would have the same minimum drajg as the circular area. 

This result was also t>btaihed by Prof. Prandtl by another method. 4 A reduction of the 
original problem of variation to the two-dimensional flow sometimes enables a survey of the 
result to be made without calculation. For instance, let a third aerofoil be added between the 
two aerofoils of a biplane having a small gap. (The gap may be about one-sixth of the 
span.) Then, in order to find the most favorable distribution of lift, the double line about 
which the flow occurs is to be replaced by three lifting lines. Now, in the region of the middle 
lifting line the velocity is small, even before this line is introduced. Therefore the discontinuity 
of the potential along the middle line is very much smaller than that along the others. Hence 
it results that the middle aerofoil of a triplane should lift less than the other two. 



« Flm communication concerning t kin Zeitachrift fur Flugtechnik und Motorl. 1814. 6. 239, in ft note by BeU. 
« Bee Techniflche Berichte der Flo*wu*mei*terftl Bd, II Heft 8. 
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PROCEDURE FOR THE CASh <F FLUIDS WITH SMALL VISCOSITY. 

The preceding results do not apply so much to the ^nJculation of the most favorable distribu- 
tion of lift as to the calculation of the least drag. For it ^opears, and the results are checked 
by calculation, that even considerable variations from the conation of most favorable distribu- 
tion of lift do not increase the drag to any great extent. Usually the minimum drag can be 
considered as the real drag of the system of aerofoils and in order to allow for the effect of 
friction of the air it is sufficient to make an addition. This addition depends chiefly upon the 
aerofoil section; it also depends, omitting the Reynolds Number, only upon the area of the 
wings and on the dynamical pressure. It is independent of the dimensions of the system of 
wings themselves. It may be useful to have a name for that part of the density of drag, inde- 
pendent of the friction of the air, which results from the theory developed in this paper. It is 
called the "induced drag/' Generally it is not the drag itself but an absolute coefficient which 
is considered. This coefficient is defined by 

Cvx-^P (52) 

where W € is the drag previously denoted by TP, q is the dynamical pressure vf.pj2, and F is 
the total area of the wings. Equation (43) can now be written 

c f « F 

where c a is the lift coefficient jr^ corresponding to c w . The greatest horizontal span 6 of the 

system of wings perpendicular to the direction of flight is arbitrarily chosen as a length char* 
acteristic of the proportions of the system, Jc is a factor characteristic of the system of aero- 
foils and has, according to the preceding, the value. 

^**^|vKj&ate.v, T ^ (54) 

y It has a special physical significance., . ^ . : , V > 

Under the same confiMto^i* single aerofoil with* a span of % times ike maximum span of a 
system of aerofoils has the same induced minimum resistance as the system. 

fcxir ftr REFINEMENT OF THE THEORY. 

The demonstrations given rest on the assumption that the velocities produced by the 
vortices are small in comparison with the velocity of flight. The next assumption, more ac- 
curate, would be that only powers higher than the first power could be neglected. 

In this case the solutions just found for lifting elements in a transverse plane can be con- 
sidered as the first step towards the calculation of more exact solutions. The following steps 

must be taken: The exact density of drag is W'^A'—-^ where v is the horizontal velocity 

produced at the lifting elements by the transverse vortices. It can be calculated exactly 
enough from the first approximation. Now, the condition of least drag is 

w-cos sin/3=*-u? 0 cos£ tt&O 

* ! and the flow of potential, according to this condition at the boundary, is to be found. Compared 
with the first approximation the density below is in general somewhat increased and the den- 
sity above is somewhat decreased. The minimum drag changes only by quantities of the 
second order. 
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If the lifting elements are distributed in three dimensions a similar refinement can easily 
be found. In this case there is to be taken into consideration a second factor which always 
comes in if the differences of the longitudinal coordinates of the lifting elements are consider- 
able. The direction of the longitudinal vortices do not agree exactly with the direction of 
flight, but they coincide with the direction of the velocity of the fluid around the aerofoil. 
They are therefore somewhat inclined downwards. A better approximation is obtained by 
projecting the lifting elements not in the direction of flight but in a direction slightly inclined 

upwards from the rear to the front. This inclination is about ~°- Except for this, the method 
of calculation remains unchanged. 



